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We address the dynamics of a drop with viscosity Xr] breaking up inside another fiuid of viscosity 

/«^ ' rj. For A = 1, a scaling theory predicts the time evolution of the drop shape near the point of snap-off 

O^ I which is in excellent agreement with experiment and previous simulations of Lister and Stone. We 

^^ ■ also investigate the A dependence of the shape and breaking rate. 

T-H ' When a fluid droplet breaks, as shown in Figure 1, a singularity develops due to the infinite curvature at the 

r^ . point of snap-ofF [Q . Near such a singularity, the axial and radial length scales become vanishingly small, allowing, 

QJ ' independent of initial conditions, a local analysis of the flow equations. Such a separation of scales implies that near 

(J^ i snap-off the profiles should be self-similar: on rescaling by the axial and radial scales the profiles near the singularity 

T^lj- should collapse onto a universal curve. [Q 

O^l ■ The character of the singularity depends on which terms in the Navier-Stokes equations are dominant at the point 

' of breakup. If the drop breaks up in vacuum, surface tension,viscous stresses, and inertia are balanced asymptotically, 

r^ ■ although the motion may pass through other transient regimes, depending on viscosity |^,^,||j|,^. In this paper, we 

^~^' investigate the situation where the viscous effects of the inner and outer fluid are included as are the pressure gradients 

' O , produced by the curvature in the surface separating them; the inertial terms are taken to be insignificant so that 

^ ' we are in the Stokes regime [0-||. Assuming that molecular scales are not reached first, this is the final asymptotic 

f^ ,; regime describing flows near snap-off for any pair of fluids even in the case of arbitrarily low viscosity. This paper 
uses experiments, simulations and theory to characterize the self similar approach to snap-off in this regime. 
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[FIG. 1 about here.; 



We consider the rupture of a fluid of viscosity Xrj surrounded by another fluid of viscosity ry. The interface between 

the two fluids has surface tension 7. At a time t* before the rupture, dimensional analysis suggests that all length 

scales have the form H{X)jri^^t* where H{X) is a function yet to be determined. Hence, if drop profiles near rupture 

are rescaled by t* , they should collapse onto a universal curve, independent of the initial conditions. However, Lister 

K^ ' and Stone M noticed that the long-ranged character of the Stokes interaction leads to logarithmic corrections in 

I , the velocity field. They simulated equations (0)-(P) below for drops having various unstable initial conditions, and 

[~^ ' demonstrated collapse if the logarithmic term was subtracted. (See also Loewenberg et al. [p|.) 

^^ , Herein, we demonstrate that this collapse also works for experiments, and construct a scaling theory to explain 

Cn ■ the profile shapes, by incorporating the nonlocal contributions into a local similarity description. Figure 2 shows the 
collapse of rescaled profiles at A = 1 for both experiments and numerical simulations using a numerical technique 
similar to that of Lister and Stone. They are superimposed on the scaling theory developed below (black line). 



^ '. [FIG. 2 about here.] 

c/2 ' The experiment used 9.5 St Glycerin dripping through 10 St PDMS. The viscosities are large enough that the 

^^, experiment is in the Stokes-flow regime even at macroscopic scales. The surface tension 7 was measured using the 
pendant drop method |10| and the viscosity was measured using calibrated Cannon-Ubbelohde viscometers. We used 
a Kodak Motion Corder Analyzer to capture ten thousand frames per second. These images were then analyzed using 
an edge-tracing program, and smoothed |11[| . 

In rescaling the experimental profiles, we shifted the origin so that the locations of the minima in the profiles lined 
^ up. Because the profiles were relatively fiat along the axial direction there was some uncertainty in the determination 

5^ , of these minima. We therefore shifted each rescaled profile in the axial, ^, direction to minimize the cumulative 
deviation in -ff(^). 

The inset of Figure 2 shows that near snap-off /i,„i„ (t) is a linear function of t* . By fitting the prefactor of this linear 
dependence, we obtain hmin — (0.031 ± 0.008)777^^t*, in excellent agreement with the result /imm = 0.0335777~^i* 
from numerical simulations 0] and the scaling theory constructed below. 

Scaling Theory Since the Stokes equation is linear, the fluid surface velocity can be expressed as an integral over 
the surface of the fiuid-fluid interface. At A = 1 the equation is [h3] 

v(^)(z,t) = -7 f K{z')3{z,z')n{z')dz', (1) 



where n is the outward normal, k is the curvature, z is the axial coordinate, and the tensor J is 
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with r the vector between the two points on the surface, I the identity matrix, and the integration is over the azimuthal 
angle 9 . Physically, equation (0) represents the response of the surface tension forcing the interface. For unequal 
viscosities A 7^ 1, eq. (0) must be amended by an additional term, which accounts for the jump in viscosity. Given 
the radial Vr and axial v^ components of the surface velocity, the interface advances according to 



dth{z,t) +Vzdzh = Vr, 

which states that the surface at a given axial position can deform by radial motion and axial advection. 
Motivated by the simulations of Lister and Stone M , we try the similarity ansatz 

h{z,t) = Vr,t*H{0, e = V^\z*/n +bllit*+^o, 



(3) 



(4) 



where z* is the axial distance from the singularity, b is a constant, and the factors of w,, = 7/77 have been inserted 
to make H and ^ dimensionless. The shift bint* in the similarity variable ^ results from the logarithmic divergence 
of the axial velocity field Q, and ^0 will be shown to be an arbitrary constant which depends on the boundary 
conditions. Since the solution near snap-off must match onto the outer profile, which varies slowly on the time-scale 
t* H{^) ^ s±£,, as ^ ^ zbcx). Here we define s_ as the slope of the shallow side of the pinch region, which by convention 
we place to the left of the minimum, and s+ as the steep slope. 

The subtle feature of this problem is the interplay of the local singularity with the nonlocal fluid response from the 
Stokes flow. The principal nonlocal effect is that the surface tension force from the cones produces a logarithmically 
diverging axial velocity field at the pinch point, C = [p[. For a local scaling theory, this singularity must be 
absorbed. We fix two points £,'_ and £,'_^_ within the linear part of the solution to the left and right of ^ = 0. Splitting 
the contributions to the velocity on the surface into a contribution from ^- < £, < £,+ and from the rest of the drop, 
and converting to similarity variables, we find 
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where V^"^' = (v^"^' — 61ni*e^)/w^ and e^ is the unit vector along the axial direction, 
axial component of the J-integral in angular brackets diverges logarithmically as t* 



-(s+(i + 4)-i + 5_(i + 



Because of the cones, the 

0. For the special choice 

s'l)~^)/4: the singularity cancels and the term in angular brackets remains finite for 



6 = 

t* — > 0. It is straightforward to extend this scaling theory to arbitrary A; in this case the amplitude of b only depends 

on A through s+, ,s_. The remaining constant in (^ depends on the detailed shape of the drop as well as on the choice 

Inserting the similarity form (0) into the equation of motion for the interface (0) gives 



fIJ-T 

dS, 



(6) 



where we have absorbed the constant advection velocity b In t* into Vz ■ 

The system (||)-(||) now has to be solved in the limit C+/^* ^ oo, (,'_/t* -^ — cx), for the interval — oo < ^ < oo 
and with boundary conditions i? ~ ^ as ^ — > cxd. Changing the constant ^o + ^ only results in a constant shift 
of the similarity function H{£^). The computation involves solving an integro-differential equation with a nonlocal 
constraint: the parameter b in (5) must be determined self consistently with the solution H{£^) according to relation 
(h). The difference between this scaling theory and others developed for fluid rupture is that here the parameters in 
the similarity equation must be determined self consistently with the solution to the similarity equation. 

We solved this system by discretizing H{£^) in an interval ^ £ [—£,inAin] and approximated all derivatives and the 
integral by second-order formulas. At ^ = — ^m,$m we demand H" = 0. Using a linear approximation for H outside 
the interval [—^in,£,in], the logarithm is subtracted explicitly. A numerical solution of the full PDE's provided an 
initial condition for Newton's iteration, which converged in a few steps. The iteration always converges to the same 
solution for any given ^o- The calculation gives Ho{X = 1) = 0.0335, s_ — —0.105, s+ = 4.81, where s+,_ are the 
asymptotic slopes at icx). These results are in good agreement with both the simulations of |§] and experiments (Fig. 



2). Although the theory has been solved only for A = 1, by continuity we expect that solutions exist for a range of A 
and that hmin obeys the law hmin = Ho{X)Vfit* . 

It is noteworthy that droplet breakup between fluids of equal viscosities is not plagued by the iterated instabilities 
found for droplet breakup in vacuum |l3]: neither experiments nor simulations observe such instabilities. The reason 
for this can be found by stability analysis of the similarity solution, following the same procedure as flj]. The result 
is that perturbations around the similarity solution can be amplified by a factor of « 150, which is much smaller 
than the corresponding amplification factor for rupture in vacuum, where it is 10^'^(i*)~^'^ [H. The reasons for the 
differences between the two problems are that (i) For A = 1 Tomitoka's formula |lq] implies that the maximum linear 
growth rate for perturbations is approximately (t*)^^. On the other hand for A = oo, the maximum growth rate is 
approximately 5(i*)~^. (ii) The A = oo problem has a time dependent amplification factor because the axial length 
scale has a different scaling with t* than the radial length scale. This implies that the closest a perturbation can be 
to the stagnation point depends on t* , even when expressed in similarity variables. 

Arbitrary A We now extend features of the above results to arbitrary A. Using Glycerin/water mixtures (1 St 
< rj < 9.5 St) and silicone oils (1 St < 7y < 600 St) we were able to cover a range of A between 0.002 and 30. The 
same procedure as used above verified self-similar data collapse in experiments with 0.02 < A < 30 [|ll|. In all these 
experiments, the conical profile associated with s_ collapsed for all analyzed profiles. This corresponds to a time 
interval of at least 0.1 seconds prior to and all the way up to the point of snap-off. This can be seen in Figure 2 for 
A = 1. In contrast, as also seen in that figure, the conical profile associated with s+ showed a time-dependent collapse 
with the region of self-similarity growing as the point of snap-off is approached. This time dependence changes as we 
vary A, slowing down as A — > 1. 

Figure 3 shows the cone angles s±, and the dimensionless breaking rate Hq as a function of A. As shown by Figs. 
3a-b, the cone angles appear to obey a power law over an extended range of A: s+ ^ ^o.2i±o.o5 g^^^j^ g_ ^ ^-o.23±o.o7^ 
(Due to measurement difficulties, the range for s_ is reduced from that for s+ UM.) Within error, the analyses 
performed on both the snap-off event near the nozzle and the snap-off event near the bulb lead to the same results. 
This agreement implies that the results are robust and independent of small variations in the surrounding flows. 
Note that our flndings are in qualitative disagreement with lubrication-type scaling arguments MM, which predict 
s± '^ A^^/^ for the slope on cither side of the minimum. On the other hand, the trends in our data are consistent 
with recent numerical simulations of the full Stokes equations by Zhang and Lister p9| . This is yet another indication 
of the breakdown of one dimensional models in describing the dynamics of two fluid rupture. 

[FIG. 3 about here.] 

A Simple Theory for Hq{X) follows by noting that (assuming the shape of the drop is slender near h„iin) the 
maximum rate that the drop can break is given by the maximum linear growth rate ri(/i,„j„) for a cylinder of uniform 
radius hmin- Namely, we have the upper bound 

By using Tomotika's formula |15[| for il{h„iin) with h„iin — Ho{X)vjjt* , this equation turns into an upper bound for 
Hf){X). This upper bound is compared with the experimental data for Hf){X) in Fig. 3c. All of the data obeys the 
bound; moreover, in the range 0.1 < A < 10 the agreement is nearly exact. Note that while most of the experimental 
data (and the upper bound equation (M) can be fit with a power law of exponent —0.53 ± 0.05, a significant trend 
with an overall negative curvature is observed in the experimental deviations. 

The agreement of the experiments with the upper bound in the range 0.1 < A < 10 is reminiscent of the marginal 
stability hypothesis, as formulated for the selection of traveling waves propagating from a stable to an unstable 
state [M. Both experiments and numerical simulations show that the breaking rate is approximately equal to the 
growth rate of linear instabilities around a cylinder of radius hmin- The upper bound in equation (M) should apply 
to all problems involving singularity formation in a system with a local instability. We have tested this upper bound 
on similarity solutions from several other examples including spherically symmetric gravitational collapse p7[ | and 
chemotactic collapse of bacteria pq] ; the upper bound is obeyed in each case, giving a reasonable estimate for the 
prefactor. Hence, this principle appears to be of rather general applicability. 

To conclude, we have (i) constructed a similarity solution for rupture at A = 1, agreeing with previous numerical 
simulations S; (ii) shown that experiments agree quantitatively with this similarity solution, both in the form of 
the profile and its time dependence; and (iii) presented a simple argument which quantitatively predicts the breaking 
rate. Experiments have also shown self-similar behavior for the range 0.02 < A < 30. There are many unresolved 
issues: Among them, there is no solid simple argument for the A-dependence of the slopes s-, s-^.. Finally, our results 



suggest that the scahng (Q) holds even in the hmit A ^ cxd, while a different set of scaling exponents is found for a 
Stokes fluid breaking up in vacuum (A = cxd) H. In addition, the profiles are asymmetric for A ^ oo, as also found 
in a recent numerical simulation EG], but are symmetric for breakup in vacuum [H, making this a singular limit. 
Preliminary experimental results of an inviscid fluid breaking up in a viscous fluid suggest that the snap-off shape is 
different from that in Stokes flow with A — > implying that this limit is singular as well. 
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FIG. 1. A drop of 9.5 St Glycerin dripping through 10 St PoIyDimethylSiloxane (PDMS) near snap-off. The nozzle diameter 
is 0.48 cm. 
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FIG. 2. The inset shows the minimum radius, hmin,{t), as a function of time for the drop shown in Fig. 1. The solid 
line is the theoretical prediction. The main figure shows the similarity function H{^) as defined by (j^. The dots are rescaled 
experimental profiles corresponding to the times indicated as arrows in the inset. The solid line is the theory, and the x's mark 
the final simulation profile. 
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FIG. 3. The asymptotic slopes s+,s-, and the rescaled minimuin radius Ho as a function of viscosity contrast A. The 
dashed Unes are a fit to the experimental data, in the right most graph the solid line is the result of our stability argument. 



